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Abstract. We investigate some examples of quantum Zeno dynamics, when a system 
undergoes very frequent (projective) measurements that ascertain whether it is within 
a given spatial region. In agreement with previously obtained results, the evolution is 
found to be unitary and the generator of the Zeno dynamics is the Hamiltonian with 
hard-wall (Dirichlet) boundary conditions. By using a new approach to this problem, 
this result is found to be valid in an arbitrary A-dimensional compact domain. We then 
propose some preliminary ideas concerning the algebra of observables in the projected 
region and finally look at the case of a projection onto a lower dimensional space: in 
such a situation the Zeno ansatz turns out to be a procedure to impose constraints. 
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1. Introduction 

Very frequent measurement can slow the time evolution of quantum mechanical systems. 
This is, in a few words, the quantum Zeno effect (QZE), by which transitions to states 
different from the initial one are gradually suppressed as the measurement frequency 
iV becomes very large PJ 12] (for a review, see pi]). There are, however, two important 
issues that deserve attention: firstly, for a general (incomplete and nonselective [I]) 
measurement, represented by a complete set of projections onto multidimensional 
subspaces (rather than a single- dimensional one, as in the usual formulation of the 
QZE, by which the measurement ascertains whether the system is still in its initial, 
pure state), the quantum system may — and indeed does — evolve away from its initial 
state, although it remains in the subspace defined by the measurement (and represented 
by a multidimensional projection operator) E]. This leads to the formation of the 
"Zeno subspaces" [7]. Secondly, if the measurement is not very frequent, the quantum 
evolution yields the so-called "inverse" or "anti" Zeno effect, by which transitions away 
from the initial state (or in general out of the relevant subspaces) are accelerated jHj- 

Both the Zeno and inverse Zeno phenomena have been experimentally observed 
during the last few years [HlUmHIllE] (but see [13] for previous analyses of experimental 
data on nuclear hadronic cascades). The first experiment was done with an oscillating 
system according to an interesting proposal by Cook jHj, and was widely debated 
[To] . In a recent beautiful set of experiments, performed by Raizen's group, first the 
initial quadratic and non-Markovian Zeno region was observed jlOj . then both the 
quantum Zeno and inverse Zeno effects were proved for bona fide unstable system 
(probability leakage out of an optical potential) [TT] . 

In this article we shall mainly analyze the first issue, investigating the features 
of the Zeno (sub)dynamics in the relevant subspace. This and related problems were 
contemplated in the seminal formulation of the QZE [2J, where it was proved that the 
dynamics is governed by a semigroup. The details of the dynamics had interesting 
and challenging mathematical aspects, that were independently investigated by other 
authors Ej. As a matter of fact, some mathematical issues are still unresolved 
nowadays. One of the most intruiguing features of the original paper [2] is that 
some delicate operator properties were postulated on physical grounds; curiously, these 
postulates are always found to be valid in concrete examples, even nontrivial ones. 

For a wide class of measurements, namely those represented by spatial projections, 
one can prove that the system evolves unitarily in a proper subspace of the total Hilbert 
space, the generator of the dynamics being the Hamiltonian with Dirichlet boundary 
conditions on the region associated with the spatial projection E]. This finding 
motivated further interesting studies on this topic [THl [THl 1201 EE] • In particular, Exner 
and Ichinose [21] analyzed this result in a rigorous framework, under the nontrivial 
(and interesting) assumption that the original Hamiltonian be lower bounded and the 
Zeno Hamiltonian densely defined in the Hilbert space. The aim of this article is to 
further elaborate on these issues. We will first explicitly work out some examples — 
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essentially the free case in two and three dimensions, with projections onto regular 
domains — and introduce a novel calculation technique, giving a constructive proof of 
the Zeno Hamiltonian. We then extend this result to a general spatial projection in N 
dimensions. 

We shall prove that the Dirichlet boundary conditions are a consequence of the 
Zeno procedure (different proofs can be given, at different levels of generality and 
mathematical rigor, see [111130121]), by exploring an interesting method of calculation, 
based on asymptotic techniques, that yields a stationary Schrodinger equation with the 
appropriate (Dirichlet) boundary conditions for its eigenfunctions. In Sec. El we set up 
the general framework and introduce notation. In Sec. El the projection domain is a 
rectangle in the plane. In Sec. [Hit is an annulus in the plane. In Sec. |S]we look at a 
spherical shell in R 3 . In Sec. HO we generalize to regular domains in M. N and in Sec. [3 
we briefly discuss the Zeno dynamics in the Heisenberg picture as well as the features 
of the algebra of observables in the projected domain. In Sec. [SI we look at a different 
case, when the system is projected onto a domain of lower dimensionality: we shall only 
look at some examples and shall not attempt to generalize. One can say that in this 
case the Zeno ansatz yields a procedure to impose a constraint. The ideas we propose 
in these last two sections are somewhat embryonic and can be considered as plans for 
future developments. In Sec. M we comment on future perspectives and applications. 



2. Zeno subdynamics 

Consider a free particle in iV-dimensions with the Hamiltonian 

2M 2M ' w v ' 

acting on if> e L 2 (R N ). Given a compact domain D C R with a nonempty interior and 
a regular boundary, consider the projection operator 

P = XD (x)= [ d N x\x)(x\, P^(x)= X d(x)^(x), (2) 

J D 

where Xd{&) is the characteristic function of the domain D, and thought of as an 
operator, along with its complement Q — 1 — P = 1 — Xd(x), decomposes the space 
L 2 (M. N ) into orthogonal subspaces. The Zeno subdynamics evolution operator is given 
by the limit 

U z (t)= lim (G(t/N)f, (3) 

where the (nonunitary) evolution 

G(t) = PU(t)P (4) 

represents a single step (projection-evolution-projection) Zeno process. 

Under rather general hypotheses the limit (J3J) can be proved to exist [111 [21] 
and yields a unitary evolution group in a proper subspace of L 2 (D). One gets 

U z (t) = Pex V (-iH z t/h), (5) 
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where the generator of the dynamics is the Zeno Hamiltonian 

Hz = ~ w (6) 

defined in the domain 

D{H Z ) = G L 2 (D) | G L 2 (D),ip(dD) = 0}, (7) 

dD being the boundary of D (hard- wall or Dirichlet boundary conditions). 

We will focus on this problem by looking for the eigenbasis {|fi)} of Uz{t) in the 
subspace PL 2 (R N ) ~ L 2 (D) such that 

(n\U z {t)\m) = lim (n\G(t/N) N \m) = (n\e- iHzt ' n \m) = 5 m . n e- lE " t/h .{8) 

N— >oo 

In order to find this basis consider an arbitrary orthonormal complete set of functions 
in L 2 {D) 

■V n (x) = (x\n) (9) 
and take the matrix elements of the single-step operator 

G m , n (t) = (m\G(t)\n) = Tr [G(t)\n)(m\] . (10) 
If the matrix elements of the single-step operator behave like 

G m>n {t) = 5 m . n (l - i^-j + R m , n {t), (11) 

where for t — > 

Rm, n (t) = 0(t), (12) 

then, under the assumption of uniform convergence of the infinite sums stemming from 
the insertion of N — 1 resolutions of the identity in (J3J), one obtains: 

G z m>n (t) = {m\U z (t)\n) 



lim Yl (t/N)G nun2 (t/N) ■ ■ ■ G nN _ un (t/N) 



n\,...,n N _ 1 



,exp -i—r- ■ (13) 



h 

The basis {|n)} is thus the eigenbasis of H z belonging to the eigenvalues E n : 

H z S> n (x) = E n ^ n (x). (14) 

Notice that when we apply U(t/N) to the relevant subspace PL 2 (M. N ), the transformed 
space need not be orthogonal anymore to QL 2 (R N ), where Q = 1 — P, and the t/N- 
dependence of the scalar product of two vectors in these two subspaces is given by 

QU(t/N)P = 0(t/N). (15) 

It has been shown that Eq. (fTTj) implies Dirichlet boundary conditions for the states 
^ n (x) (the "Zeno eigenbasis") (BJ. The proof, based on asymptotic techniques, yields 
the propagator in an appropriately chosen basis of eigenf unctions. In the following 
sections we shall introduce a novel approach: by using asymptotic analysis and the 



Zeno dynamics and constraints 



5 



path integral representation of the matrix element (|1(J|). we will obtain a stationary 
Schrodinger equation and a set of boundary conditions for its eigenfunctions. This 
will enable us to define the induced Zeno Hamiltonian Hz and its spectrum. The 
advantage of the present approach, as compared to the previous one 0, lies in the 
fact that one can derive a Schrodinger equation with Dirichlet boundary conditions 
in the projected (Zeno) subspace. Moreover, by examining some examples of multiple 
connected domains, we will show how the Zeno dynamics induces constraints that inherit 
the topological properties of the parent space. 



3. Rectangle 

We start off with one of the simplest examples and introduce the procedure. Consider 
a rectangle in the plane, D = [0, a] x [0, b] C M 2 . In this case the projection (j2j) reads 

pa rb 

p = xm( x ) xio,b](y) = / dx / d v\ x y)( x y\ ( 16 ) 

Jo Jo 

and the Hamiltonian is 

H - P W-m^- < 17 > 

The Zeno Hamiltonian, engendering the Zeno subdynamics, is formally given by ©-© 
and represents a free particle in the box D = [0, a] x [0,6] with Dirichlet boundary 
conditions 

H Z = -^(dl + d 2 y ), (18) 
^(0, y) = V(a, y) = 0, $(x, 0) = ^{x, b) = 0. (19) 

The eigenfunctions and eigenvalues are well known 

*nm(a?,2/) = sin (^-^) y| sin i^v) , ( 20 ) 
h 2 ir 2 ( n 2 m 2 



Enm = ^w{^ + ^)- (21) 

Let us look in detail at the derivation of the Zeno subdynamics (I18|) - (j21|) in this 

particular case. As explained in Sec. |2l the eigenbasis of the Zeno Hamiltonian Hz 
in L 2 {D), 

^nm(x,y) = (x,y\nm) (22) 



Gn'm',nm(t) ~ ^n'n^m'm [ 1 ~ i T~ ) + o(t), (23) 



must satisfy condition (fTTj) : 

A. , ,Y ,. 11 -i- 

h 

where 

Gn'm' ,nm 

(t) = (n'm'\G(t)\nm). (24) 
are the matrix elements (jl(Jj) of the single-step evolution operator. 
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This can be proved by direct inspection: one gets 



dx l dy l dx 'L dy i^L) 



x^* n , m ,(x,y)^ nm (x',y') (25) 
and by substituting £ = x' — x and r\ = y' — y 

fb ra—x rb—y / 



ra rb ra-x rb-y f M \ 

G nlm , nm (t) = J q dxJ Q dy J ^ d{ J drj 



g 2ht 



-y 

x^* nlml (x,y)^ nm (x + ^y + r]). (26) 

With the natural choice ty nm (x,y) = ip n (x)(j) m (y) this yields the product of two 
quantities 

G n 'm',nm(t) G n ' n (t)G m ' m (t) 

a pa—x / j^j \ 1/2 ^ 2 

dx I d t\7TZ7U) e i ^ip^{x)ip n {x + ^) 



2niht 



o J -x 

b fb-y / M ^ 



xj o dy'J dv[^jz) e^cl>* m ,(y)My + v) (27) 

and accordingly E nm = E n + E m . Consider the first quantity G n > n and the integral over 
£. In the small-t limit there are contributions from the boundary points £ = —x and 
£ = a — x and from the stationary point £ = 

r 

G n ' n = / dxip n , (x) [bound + stat] , (28) 
Jo 



where 



bound = 7T-T+ -J7^^ X + ^ e 2ht 
2m fit J 



£=a— X 



+ O (t 3/2 ) 

I fcl / p iM{x-a) 2 /2M AMx 2 jlU \ 



(29) 



while (A = M/2ht) 



stat = / de ( v>„(ao + <(x)e + ^:(x)£ 2 + o(e)) J ^-/ xe 



th 

= Ux)+^'n{x) + 0(f). (30) 
In order to obtain one must require that (remember that E nm = E n + E m ) 

bound = 0(t 3 / 2 ) and - ^<(x) = ^ty»0«0, (31) 
which translates into 

" 2M^ n(x) = En ^ X ^ With ^n(0)=^„(a) = 0, (32) 
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so that for G n > n one obtains 



G n , n (t) = (i_^)5 n , n + 0(t 3 / 2 ) (33) 



and analogously for G m ' m , so that 



E t E t s 

G n 'm',nm(t) = I 1 — « — ^~ ^n'n^m'm + 0(£ 3/2 ), (34) 



h h 

which has exactly the form (J23J). By Eq. (J32j) and its analog for 4> m (y), the eigenfunctions 
*nm(i,2/) = ^n{x)<p m {y) of #z satisfy 
ft 2 

-^(^ + a 2 )^ n , m (x,y) = (£ n + £ m )^ n , m (x,2/), (35) 

with Dirichlet boundary conditions. They are therefore given by (|20p. The Zeno 
Hamiltonian is therefore (|T8|) - (|T9|) . 

This derivation, although it yields the desired (and correct) result, is not rigorous. 
In particular, it does not tackle the delicate problem of understanding the convergence 
properties of the asymptotic expansion at the intersection of the (x, y) and (a/, y') 
boundaries in Eq. [this is apparent if one looks at the denominators of the far 
right hand side of Eq. (|29j) ]. A similar approach will be adopted in the next sections. A 
more rigorous proof can be given, but will not be presented here. 

4. Annulus 

Consider now a circular annulus (or ring) of width 5r = r 2 — T\ on the plane, defining 
the domain D = {(x,y) \ r\ < x 2 + y 2 < r 2 }. The projection on D reads 

r pr2 r'i-K 

P = X[n,r 2 ]( r ) — / dxdy\xy) {xy\ = / drr / dQ\rO)(r6\. (36) 
Jd Jn Jo 

r 2 — t\ = 5r > 0. (37) 

The Zeno Hamiltonian, engendering the Zeno subdynamics, is given by © and 
represents a free particle on D with Dirichlet boundary condition 

i/.(r 1 ,e)=V'(r 2 ,9) = 0. (39) 

As is well known, by writing the eigenfunctions of Hz as ^ n i(r, 9) = ip n i(r)(fii(6), the 
angular functions are given by 

(j>i{9) = -=exp(iZ0), with I = 0, ±1, ±2, . . . , (40) 
v 2ir 



while the radial part of the eigenvalue equation reads 

h 2 l ft 2 / 2 

-d r (rd r )ilj nl (r) + TTTT^^mir) = E nt ij; n i(r), (41) 



2Mr ~<v~< / r™v, ■ 2Mr 2 

«?~i) = V^K> = (42) 
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and can be solved in terms of Bessel functions. 

Let us look in detail at the derivation of the Zeno subdynamics f)38p -()42j ) in this 
particular case. As explained in Sec. El the eigenbasis of the Zeno Hamiltonian Hz in 

*ni(r,9) = (r,9\nl) (43) 
must satisfy condition (fTTf . that is 



GW,»i(t) = Sn'Jin U - i^f) + o(t), (44) 

where 

G nll , >nl {t) = (n'l'\G(t)\nl) (45) 

are the matrix elements (fTU|) of the single-step evolution operator. 
By writing ^ n i(r,9) = ip n i(r)(pi(9), we get 

T2 /T2 /*27T /»27T 



/T2 /*7*2 /*^7T /*^7r 

Gf„/,/ ini (t) = / rrfr / r'dr' / d0 / d0' 

Jri Jrx JO JO 



x MrMWrtvWti'W J e'^S", (46) 

where d is the distance between the points (r, 0) and (r', #') 

d 2 = r' 2 + r 2 - 2r'r cos(9' - 0) = (r - r) 2 + 2r'r(l - cos(0' - 0)). (47) 

Let us look first at the 6 integrals. Changing again to rj = 9 — 9' and dropping the prime 
one gets 

r 2 r 2 ( M \ 1/2 -Mtr'-r) 2 

G n n>,nl{t) = J rdrj r'dr> nZ (r )<,,(/) [^-^ t ) e*^^ 

(48) 

Consider the integral over rj (at fixed r' and r). In the limit t — > the boundary 
contribution reads (z = Mr'r/ht) 

bound = j_ . 7\ ^ft-cosffl) ^ (27r) _ 0/(o)] + 0(f 3/ 2) (49) 
V r'r V2niz sin 

In order that 0(\/i) vanishes and (jUj) is satisfied, one must require the periodicity 

0,(0) =<M2vr). (50) 

The difference with the preceding case is given by the periodicity of the Green function. 

However, we now have two stationary points in the r\ integral. One is rj = and 
the other is rj = n for 9 < 7r, or rj = —n for 9 > n. These represent the minimum 
and maximum of the distance between two points, one fixed on the circle r' = const, 
and the other one located on the circle r = const, at an angle 77. One should get (at 
least) two points of stationary phase each time one constrains the dynamics on a closed 
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(iper) surface. Both contributions must be taken into account. The only difference with 
the previous case is that one must consider also n 4 terms arising from the cosine in the 
integral 

1 



statr 



st at -J-7T- 



dr\ 



2m 



^izrf /2 



2iz 



( j> l (e±7r) + - ( f)' l '(e±7r)ri 2 

-izif/2 



(51) 



(52) 



Notice that static has a phase 2z = 2mr'r/ht. This phase changes the term 
m(r' — r) 2 /2ht of the r',r integrals into a term m(r' + r) 2 /2ht. This factor has no 
more stationary points in the r', r integrals, so that its contribution can be neglected 
(in the t — > limit). In turn, also the contribution from stat-t,,- can be neglected. On 
the other hand, the stato contribution is 

1 



statr 



+ 0(t 2 



(53) 



IMr'r 1 v ' 8Mr'r 

and following the same reasoning as in Sec. 121 (rectangle on the plane) one obtains a 
differential equation for the eigenfunctions 

-#(0) = a,&(0), 0(O) = 0(2vr) (54) 

which yields ct\ = I 2 , whence 



2- 



d8(J)y(8) stat = 5i>r 
Therefore the integral over r', r reads 



ti- 



ll 



2Mr'r 



I 1 



0(t 2 



(55) 



r-2 



rdr 



r-2 



r'dr' 



M 
2iriht 



■ 111 {r —r) 

iPn>i>{r)ipni{r)e l m 8 n 



x 



h 




1 - it 


('-91 


2Mr'r 





(56) 



By inserting £ = r — r' and dropping the prime on r' we get 

T2 fTi—r 



n 



X 



y/rdr 
1 - it 



r\ —r 



M 
2nifit 



h 



2Mr(r + f ) 



i 1 -- 



(57) 



By the same reasoning as before one obtains a differential equation and the Dirichlet 
boundary conditions for the functions A n i(r) = y/ri/i n i(r): 



2M 



h 2 



2Mr 2 



l 2 --jA nl (r) = E nl A nl (r), (58) 

A nl {n) = A nl (r 2 ) = 0. (59) 

In terms of the radial functions ip n i, Eq. ()58|) becomes just Eq. (}4*T]) . whence the Zeno 
Hamiltonian is given by (|HS |) -(|H9 ^) . 
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It is interesting to notice that in this case of multiple connectedness the Zeno 
dynamics yields no Aharonov-Bohm topological phases. In a few words, one might say 
that the projected dynamics on the annulus "inherits" the topological properties of the 
initial space M 2 , and in particular the single valuedness of the wave function. The spatial 
projections do not introduce any additional "twist" into the system, that could induce 
a phase. 

Two additional quick comments: first, the r\ — > limit yields a circle; however, it 
does not yield the Zeno dynamics on the domain D = {(x, y)\x 2 + y 2 < r|}, because 
of the spurious condition i[) n i(0) = 0, excluding s-wave eigenf unctions. This seemingly 
trivial remark clarifies that taking a limit of the projected domain does not necessarily 
yield the right Zeno dynamics. Second, the circular ring sector {(r, 8)\r 1 < r < r 2 ,6i < 
6 < 62} can be easily computed and yields the right dynamics and eigenf unctions (Bessel 
functions I^r), \i <E K) |22J - 



5. Spherical shell 

Let us now consider an example in M. 3 . We first observe that the parallelepiped 
can be easily dealt with by extending the techniques of Sec. El We therefore look 
at a more interesting situation. Consider a spherical shell in M. 3 and a domain 
D = {(x, y, z) I r\ < x 2 + y 2 + z 2 < rf }. The projection on D reads 

P = X[ri,r 2 ]( r ) = / dxdydz\xyz)(xyz\ 

PT2 /*7T rill 

= / r 2 dr / sin0d0 / d<p\r6<p)(r6<p\. (60) 

Jr 1 JO JO 

The Zeno Hamiltonian, engendering the Zeno subdynamics, is given by © and 
represents a free particle in the spherical shell D with Dirichlet boundary condition 

Hz= -^d 2 x + d 2 y + dl) (61) 

h 2 f„, 2 



d r (r 2 d r ) + —d e (sm6d e ) + —^di 



2Mr 2 V sinO ° K °' sin 2 9 

^(r 1 ,0,0) = ^(r 2 ,0,0) = O. (62) 

As is well known, by writing the eigenfunctions of Hz as ^ n i m (r, 8, 0) = 
Rni( r )Yi m (9)Q m ((f)) , the radial part of the eigenvalue equation reads 



21[ r2 d r (r 2 dr)Rm(r) + —^-LB^ir) = E nl R nl (r), (63) 

Rni(n) = R nl (r 2 ) = (64) 

and can be solved in terms of spherical Bessel functions. 

Let us see how one can obtain Hz in this case. The first steps of the derivation are 
the same as before. By rewriting the distance d(r'6'<j)', r6<p) as 

d 2 = (r' - r) 2 + 2rV(l - cos(fl' - d)) + 2r' rsin 9' sin 9(1 — cos(0' — 0))(65) 
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it is apparent that the integrals must be performed in the order <fi — > 9 — > r and that 
only those stationary points that do not give an additional phase contribute to the final 
result. 

As eigenfunctions we choose the orthogonal set 

VnlmW) = Rnl(r)Y lm (0)® m (<t>)- (66) 

The transition element is 

G n 'i' m ',nml(t) = / r /2 dr' / r 2 drR* n , v (r')R nl (r) 



M - M(r'-r) 2 1 

XA/ 7T~I7 e 2M -T G l'm',lm, (67) 

lixifit rr 



rn fir 

G Vm . M = / smew / smedeY* m ,(e')Y lm {e) 

Jo Jo 



x ^ Mr'r ci M r , r{1 _ cos{e ,_ e)) 
2iriht 

1 JjL/ f 2n J ! .\ f r> mu f)' r f) 



Vsin 0' sin J J V 2iriht 

x e^ r ' rsine ' sine ( 1 - cos ^'-^)$^(0 , )$ m (0) (68) 

The 4>', 4> integral is immediately computed as in the case of the annulus, Sec. HI $ m 
must therefore satisfy the differential equation 

- ®l = a m $ m , $ m (0) = $ m (27r), (69) 

so that a m = m 2 . Then Gi' m 'i m becomes 



G Vm , M = [ vWflW / Vsh[6d6Y; m ,(6')Y lm (0) 
Jo Jo 



Mr'r 
2iriht 



x e*^ 1 — ( 1 - i — ; ^ . fl (m 2 - 1/4)) <W 
\ zMr sm 0'r sin 6 J 

(70) 

The integral over 0', can be computed in a standard way (do not forget the £ 4 term 

in the cosine series) and this in turn requires that the function A[ m = \J sin BY\ m must 

satisfy the differential equation 

I m 2 -l/4 

4 sin 2 

or, equivalently, 

| Q / \ 7J1 2 

sin0— Y lm r^r7{Yi m = -ai m Y lm , (72) 



A" A — ' A A f71 "\ 



sin 00 V <90y sin 2 0' 
with Y/ m (0) = Yj m (7r). This is the standard equation for spherical harmonics. It is 
known that ai m = 1(1 + 1) irrespectively of the value of m. We obtain 



]\/[ M(r f — r)^ 

G n > Vm >nim= J r'dr'J rdr R* n , v (r')R nl {r) ^ i^-f 1 2Rt 

x f 1 " ^m^*) 5l ' l5m ' m - (73) 
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h 2 k 2 

Finally, the differential equation for A n i = rR n i reads (here E n i = 2R ™ 1 , which is 
independent of m) 

- A" nl + l ^±^A nl = k 2 nl A nl , A nl ( n ) = A nl (r 2 ) = 0, (74) 



r 



or, equivalently, in terms of R n i, Eq. The Zeno Hamiltonian is therefore given by 

6. The general case 

By looking at the preceding examples one might think that the method introduced 
in this article is parochial and works only, for example, when the domain, besides 
being sufficiently regular, is also endowed with particular symmetries (regular polygons, 
circles, spheres and so on), that enable one to introduce coordinates with a range of 
integration that can be reduced to a product of intervals. In turn, this might appear as 
an implicit condition of separability, e.g. in the case of the three-dimensional Schrodinger 
equation [22]. On the contrary, as will be shown in this section, the method we propose 
is of general applicability. 

Consider again the Hamiltonian and the projection ({2~|). D C M. N being a 
compact domain with nonempty interior and a regular boundary. The iV-dimensional 
propagator (fTUj) reads 

G m , n (t) = (m\G(t)\n) 

f „ f M ( M \ N/2 . M(x-„) 2 

= L dX L dy {^t) ^^MV) (75) 
and by substituting £ = y — x one gets 

r r f m \ N/2 Mt 2 

G m>n (t) = J d d N x ** m (x) J d d N i ( — J e 4 w* n (a. + £) 

= / d N x *™(cc) [bound + stat], (76) 
Jd 

where 

D-x = {y\x + y e D}. (77) 

Let us evaluate separately the two contributions in the small-t limit. In order to compute 
the boundary term, we first observe that 

^ = (78) 



and then integrate by parts (A = M/2ht) 

N ' 2 £ ■ Ve^ 2 



bound = J d N £l—j #„(a; + £) 



N/2 



d N iv 

D 



2i\e 
2i\e 
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D 



2 



N/2 



(—) 

\2mht ) 



N/2 



N _ lo v n (y)(v -x)-u e ^(-y) 



, dD (y - xf iM/ht 

(1 + 0(t))] , (79) 

u being the unit vector perpendicular to the boundary. In the second equality, Eq. (f7S|) 
was used again in order to obtain a higher-order volume integral with the same structure 
as the initial one. The stationary contribution is obtained, as usual, by expanding the 
integrand function around x 

stat= !d N ie^ 2 
\27TihtJ J s 

x (y n (x) + W n (x) ■ £ + ^didjVnWtej + 0(|^| 3 )^) . (80) 

Observe that the contributions of the linear and quadratic (with i ^ j) terms in the 
integral vanish due to symmetry and one is left with 

th 

st a t = V n (x) + i—AV n (x) + 0(t 2 ). (81) 

In order to obtain (jll|) - ()12)) from (|7fjjl one must require that the leading contribution in 
the boundary term (|79|) vanishes and 

' l Am n (x) = ^* n (aj), (82) 



namely 



2M y ' h 

1,2 A* w (aj) = E n * n (x), with V n (dD) = 0. (83) 



2M 

Notice that the Schrodinger equation is obtained from the stationary contribution to the 
asymptotic expansion, while the Dirichlet boundary conditions are a consequence of the 
requirement that the boundary term ()79jl vanish at the lowest order in the expansion. 

Let us briefly comment on the features of the method introduced. As already 
emphasized at the end of Section 01 this analysis, although not entirely rigorous, yields 
the correct result. We derived the desired properties of the propagator by requiring 
at the same time the validity of the Schrodinger equation and the Dirichlet (hard- 
wall) boundary conditions for the eigenbasis of the (Zeno) Hamiltonian. We should 
emphasize, however, that the boundary and stationary terms are being dealt with 
separately. In fact, we did not consider the contribution of those boundary points 
that are also stationary points. Such points belong to the intersection of the boundaries 
of the two domains D in (|7Hj) . namely x = y G dD, and should be analyzed with great 
care. A more rigorous treatment can be given, in which the contribution of the integral 
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(|76|) is uniformly estimated: this analysis requires a different evaluation of the boundary 
terms and will be presented elsewhere. 

The introduction of a potential [H] is not difficult to deal with if the detailed 
features of the convergence (fTT H) -(fTT ]) are not worked out. Much additional care is 
required at a deeper mathematical level, when the self adjointness of the Hamiltonian 
is called into question and must be explicitly proved. If additional rigorous results 
[21 EH UH I2J] are taken into account and, by an educated guess, extended to the case of 
a sufficiently regular potential, one is tempted to assume that the procedure sketched 
above is valid in general and the Zeno dynamics governed by a self-adjoint generator 
(and a unitary group). The situation may clearly become more complicated when the 
potential is singular and/or the projected spatial region (or its boundary) lacks the 
required regularity. 



7. Zeno dynamics in Heisenberg picture 

In this section we would like to consider the Zeno dynamics in the framework of the 
Heisenberg picture. The following discussion must be considered preliminary: additional 
details and a broader picture will be given in a forthcoming paper. An interesting and 
natural question concerns the algebra of observables after the projection. This is not a 
simple problem. One can assume that to a given observable O before the Zeno projection 
procedure there corresponds the observable POP in the projected space: 

O => POP. (84) 

For example, if one starts in K. and projects over a finite interval P = Xi{ x ) {I being an 
interval of R), the momentum and position operators become 

{id T for x E I , „ 

(J otherwise 

x^PxP=\l { °l XeI . (86) 
[ otherwise 

In this respect it is easy to see that the correspondence (|8*4*jl is not an algebra 
homomorphism. However, if we redefine a new associative product in the algebra of 
operators, by setting 

A* B = APB, (87) 

with this new product the previous correspondence (J54*|) becomes an algebra 
homomorphism |24j . Notice also that the new (projected) algebra acquires a unity 
operator P. 

Notice that in general the evolution will not be an automorphism of the new 
product. However, it will respect the product to order 0(t/N) and induce, in the limit, 
a Zeno dynamics on the projected algebra, i.e. on the image of the projection. The 
evolution will be trivially an automorphism when it commutes with P and is therefore 
compatible with the new product without any approximation. For instance, this would 
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be the case if we take as Hamiltonian the square of the angular momentum in the case 
of the annulus (Sec. HJ). 

In general one has to modify the associative product in such a way that the 
"deviation" of U(t/N) from being an automorphism is of order o(t/N), so that in the 
limit Uz{t) will be an automorphism of the new associative product adapted to the 
constraint. In other words, the sequence of evolution operators 

V N (t) = G(t/N) N = {PU(t/N)P) N , (88) 

yielding the Zeno limit is mirrored at the level of the algebra by the following 
sequence of deformed associative products 

A * N B = AP N B, (89) 

where P/v is a smooth positive operator with < P/v < 1 and PnP = PPn = P- 
For any N, P N forms together with Q N = 1 — P N a positive operator valued measure, 
yielding a resolution of the identity, i.e. Pn+Qn = 1, which approximates the orthogonal 
resolution P + Q — 1, in the sense that 

P Nl p = P^ + 0(1/ N), V</> e L 2 (R N ). (90) 

For any N the evolution V/v(£) is an automorphism of the product *n and in the limit 
N — > oo we get the desired result (|S7jl . 

Observe that, for unbounded operators, (}8l|) does not necessarily yield self-adjoint 
operators: for example, after the Zeno procedure, the momentum p would act on 
functions that vanish on the boundary of / and would have deficiencies (1, 1), see [5]. On 
the other hand the Zeno Hamiltonian @ is self adjoint. However, it would be arbitrary 
to require a similar property for every observable in the algebra. We shall analyze this 
issue in greater detail in a future article. In general, the lack of self-adjointness of the 
operators representing the "observables" of the system in the projected subspace might 
be related to the incompleteness of the corresponding classical field |25] 15]. 

8. Projections onto lower dimensional regions: constraints 

In all the situations considered so far, the projected domain always has the same 
dimensionality of the original space (M n ). [Remember that, after Eq. (P), we required 
the projected domain D to have a nonempty interior.] However, it is interesting to 
ask what would happen if one would project onto a domain D' of lower dimensionality. 
This is clearly a more delicate problem, as one necessarily has to face the presence of 
divergences. It goes without saying that these divergences must be ascribed to the lower 
dimensionality of the projected domain and not directly to the convergence features of 
the Zeno propagator [2^]. Our problem is to understand how these divergences can be 
cured. One way to tackle this problem is to start from a projection onto a domain 
D C M n and then take the limit D — > D' C M n_1 , with a Hilbert space (Zeno subspace) 
L\D'). 
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The content of this section is preliminary. We shall only sketch the main ideas and 
postpone a thorough analysis to a forthcoming paper, in which the physical meaning of 
the divergences will be spelled out in greater details. 

8.1. From the rectangle to the interval 

Let us first look at the case of the rectangle, investigated in Sec. |3J and let b — > 0. We 
first notice that in order to get a sensible result one must first perform the Zeno limit 
N —>■ oo and then let b — > 0. In particular one must require 

2Mb 2 

St = t/N < h/E m = -j—^, (91) 

which has an appealing physical meaning: the time during which the particle evolves 
freely between two projections must be small enough so that the particle remains well 
within the rectangle of width b. In practice, one must first set m < m*, for some m*, in 
order to obtain a sensible result and finally let m* become arbitrarily large. The order 
in which the two limits (N — > oo and b — > 0) are taken is relevant and significant from 
a physical perspective: one must first make sure that the wave function does not "leak" 
out of the projected (Zeno) region and then let this region "shrink" into a domain of 
lower dimensionality. 

However, even if one follows the correct procedure (i.e., first iV — > oo and then 
b ^ 0) one still gets divergences in the phases, since 

„ h 2 7r 2 m 2 „ , . 

E m = — -»• oo for b -> 0. 92 

2Mb 2 K ' 

Notice also that since the energy differences between different m states diverge, a 

superselection rule arises. Different subspaces, labeled by different values of the 

quantum number m, remain separated (at least for low-energy processes with energies 

E h 2 /Mb 2 ). This is also physically revealing. 

On the basis of the above insights, we therefore propose to perform the limit 6^0 

by choosing a particular eigenstate 4> m (y) and considering the reduced evolution 

U z (t) = e iEmt/h (m\U z (t)\m} } (93) 

which operates only on the x degree of freedom. Physically, this corresponds to the 
choice of a particular value of the superselection charge. Thus, the reduced propagator 
reads 

G(x',x;t) = (x'\U z (t)\x) = e iEmt/h (m;x'\U z (t)\m;x) 



E 



e - iB -" B V„(i')«(i), (94) 



where ip n are the eigenfunctions of the Dirichlet problem (JH2J), and one gets 

H z = —ojf, with Dirichlet b.c. (95) 

This is just the free particle on the interval [0,a], as expected. Not only can the 
divergence be cured, it also yields the desired result. 
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The framework explained in this particular example works also in more complicated 
circumstances. In particular, it is important to understand in which order the two limits 
must be computed: first one makes sure that the Zeno mechanism works efficaciously, 
then takes the desired limit on the domain. We consider here two other simple situations. 



8.2. From the annulus to the circle 

Let us now look at the annulus, investigated in Sec. |H We would like to recover 
the evolution of a particle on a circle by considering the 5r —>■ limit, while keeping 
fi — f~2 — R constant. Once again, as in Sec. 18. 1[ we have to face some divergences. By 
taking the limit one finds the approximate eigenfunctions of Eq. ()41j) 



2 /727T 

Mr) ~?p n (r) = \/ — sin [—(r - R)) (96) 



and the energies 



E n ^E n + El = + # (V - i ). (07) 



2M 5r 2 2MR 2 \ 4, 

Again one finds a diverging energy which must be regularized. However a second (finite) 
term appears (—h 2 /8MR 2 ) j27j which is not present in the usual circle quantization. 
We notice that different quantization methods yield different results [28]. 
The reduced propagator on the remaining degree of freedom 9 is just 

2 2 2 

G{9',9;t) = (n,9'\U z (t)\n,9) = J^^^ii ')^), (98) 

i 

which is what one expected. 
8.3. From the shell to the sphere 

Finally, we reconsider the spherical shell of Sec. El and take the limit 5r — > 0, while 
keeping r\ = r 2 = R constant, like in Sec. 18.21 This yields the energies 

h 2 n 2 7T 2 H 2 l(l + 1) 
nl ~ 2^ + 2MR 2 { ' 

and following the same regularization procedure as before we find 

G(6\ 0', 9, <j); t) = e l ^\n; 9', <f>'\U z (t)\n\ 9, <j>) 

I in 

which is the usual propagator on a sphere of radius R, whose Hamiltonian is 

i =4=-2if(^* )t i4 < ioi) 
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9. Concluding remarks on potential applications 

We have investigated the quantum Zeno dynamics, when a free system undergoes 
frequent measurements that ascertain whether it is within a sufficiently regular spatial 
region. The evolution in the projected (Zeno) subspace is unitary and the generator of 
the Zeno dynamics is the Hamiltonian with hard-wall (Dirichlet) boundary conditions. 
In general, this procedure leads to the formation of the "Zeno subspaces" [7], on whose 
boundaries the wave function must vanish (Dirichlet): this is the ultimate reason for 
the absence of amplitude (and probability) leakage between "adjacent" subspaces. 

Quantum computation [29J is one of the most promising fields of potential 
application of the QZE. Interactions with the environment deteriorate the purity of 
quantum states and represent a very serious obstacle against the preservation of quantum 
superpositions and entanglement over long periods of time. It is therefore of great 
interest to endeavor to understand whether decoherence can be controlled and eventually 
halted [SO]: in this context, novel techniques hinging upon the quantum Zeno effect are 
of interest. Besides the use of quantum error correcting codes [31], the engineering of 
"decoherence-free" subspaces is also recently being considered and widely investigated 
[32] • Some mechanisms are actually being proposed, based on the so-called "bang-bang" 
evolutions and their generalization, quantum dynamical decoupling [33] • Although 
"bang-bang" techniques in classical control theory are know to engineers since long ago 
[31], their introduction as a quantum control and their unification with the basic ideas 
underlying the quantum Zeno effect are quite recent [35]- In particular, the decoherence- 
free subspaces are the dynamically generated quantum Zeno subspaces [7j within which 
the dynamics is far from being trivial, as has been discussed in this article. It is 
also worth noticing that the range of applicability of these ideas is wide, as frequent 
interruptions and continuous coupling J3H] can yield similar physical effects. This is 
not entirely surprising [31], in view of Wigner's notion of "spectral decomposition" 
[38] . However, when one considers applications of the Zeno dynamics in the context 
of decoherence-free subspaces, one must remember that if the measurement is not very 
frequent, the quantum evolution yields the so-called "inverse" or "anti" Zeno effect, by 
which transitions out of the decoherence-free subspace is accelerated |B]. 

In conclusion, it is interesting to notice that an issue that was considered as purely 
academic until a few years ago, has been first experimentally demonstrated and is now 
being considered as a possible strategy to combat decoherence, with interesting spinoff's 
and very practical applications. 
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